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ABSTRACT. The power spectrum of mass density fluctuations is estimated 
from the Mark III and the SFI catalogs of peculiar velocities by apply- 
ing a maximum likelihood analysis, using parametric models for the power 
spectrum and for the errors. Generalized CDM models with and without 
COBE normalization are used. The applications to the two different data 
sets give consistent results. The general result is a relatively high ampli- 
tude of the power spectrum, e.g, at k = 0.1/iMpc - 1 we find P(k)Q 12 = 
(4.5 ± 2.0) X 10 3 (/i^Mpc) 3 , corresponding to crgH ' 6 = 0.85 ± 0.2. Model- 
dependent constraints on combinations of cosmological parameters are ob- 
tained for families of COBE-normalized CDM models. These can roughly be 
approximated by QhfjQ^n" = 0.6 ± 0.2, where pi = 1.3 and v = 3.7, 2.0 
for flat ACDM models with and without tensor fluctuations respectively. For 
open CDM, without tensor fluctuations, the powers are p = 0.9 and v = 1.4. 
The quoted error-bars reflect the 90% formal likelihood uncertainty for each 
model and the variance among different models and between catalogs. This 
is a brief review of a collaborative project (for more details, see Zaroubi etal. 
1997, Freudling etal. 1998). Preliminary constraints in the Q — S1a plane are 
presented as well. 



1 Introduction 

In the standard picture of cos- 
mology, structure evolved from small 
density fluctuations that grew by 
gravitational instability. These ini- 
tial fluctuations arc assumed to have 
a Gaussian distribution character- 
ized by the power spectrum (PS). On 
large scales, the fluctuations are lin- 
ear even at late times and still gov- 
erned by the initial PS. The PS is 
thus a useful statistic for large-scale 
structure, providing constraints on 
cosmology and theories of structure 
formation. 

The galaxy PS has been estimated 
in recent years from several redshift 
surveys (see reviews by Strauss & 
Willick 1995; Strauss 1998). Alterna- 
tively, one can estimate the PS us- 
ing measurements of peculiar veloci- 
ties, which are directly related to the 
mass density fluctuations. Velocities 
are also sensitive to larger scales and 
thus subject to weaker non-linear ef- 
fects. In this work, we develop and 
apply a likelihood analysis (first pro- 
posed by Kaiser 1988) in order to 
estimate the mass PS from peculiar 



velocity catalogs. The method, act- 
ing on the "raw" peculiar velocities 
without additional processing, uti- 
lizes much of the information con- 
tent of the data. It takes into account 
properly the measurement errors and 
the finite discrete sampling. The sim- 
plifying assumptions made are that 
the velocities follow a Gaussian dis- 
tribution and that their correlations 
can be derived from the density PS 
using linear theory. 

Two catalogs are used for this 
purpose. One is the Mark III cata- 
log of peculiar velocities, a compi- 
lation of several data sets, consist- 
ing of roughly 3000 spiral and ellip- 
tical galaxies within a volume of ~ 
80h _1 Mpc around the local group, 
grouped into ~ 1200 objects (Willick 
etal. 1995, 1996 1997). The other is 
the recently completed SFI catalog, 
a homogeneously selected sample of 
~ 1300 spiral field galaxies, designed 
to minimize effects of combining dis- 
parate data sets (Haynes et al. 1998; 
Wegner et al. 1998). In both catalogs, 
the typical relative distance errors of 
individual galaxies are 15 — 20%, and 
both data sets are carefully corrected 
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for the various systematic biases. It 
is interesting to compare the results 
of the two catalogs, especially in view 
of apparent discrepancies in the ap- 
pearance of the velocity fields (e.g, 
da Costa etal. 1996, 1998). 

2 Method 

Given a data set d, the goal is to 
estimate the most likely model m. 
Invoking a Bayesian approach (and 
assuming a uniform prior), this can 
be turned to maximizing the like- 
lihood function C = V(d\m), the 
probability of the data given the 
model, as a function of the model pa- 
rameters. Under the assumption that 
both the underlying velocities and 
the observational errors are Gaussian 
random fields, the likelihood func- 
tion can be written as 

C = [(2 7 r) jv det(i?)]- 1 / 2 

xexp X)*^ 1 ^ J ' 

where {di}?L 1 is the set of observed 
peculiar velocities and R is their cor- 
relation matrix. R involves the theo- 
retical correlation, calculated in lin- 
ear theory for each assumed cosmo- 
logical model, and the estimated co- 
variance of the errors. 

The likelihood analysis is per- 
formed by choosing some paramet- 
ric functional form for the PS. Going 
over the parameter space and calcu- 
lating the likelihoods for the differ- 
ent models, one finds the PS param- 
eters for which the maximum like- 
lihood is obtained. Confidence lev- 
els are estimated by approximating 
— 21n£ as a \ 2 distribution with re- 
spect to the model parameters. Note 
that this method, based on pecu- 
liar velocities, essentially measures 
f(Q) 2 P(k) and not the mass den- 
sity PS by itself. We extensively test 
the method using realistic mock cat- 
alogs, designed to mimic in detail the 
real catalogs (Kolatt etal. 1996; El- 
dar etal. 1998). 

We use several models for the PS. 
One of these is the so-called T model, 
where we vary the amplitude and the 
shape-parameter T. The main anal- 
ysis is done with a suit of general- 
ized CDM models, normalized by the 



COBE 4-year data. These include 
open models, fiat models with a cos- 
mological constant and tilted mod- 
els with or without a tensor compo- 
nent. The free parameters are then 
the mass-density parameter f2, the 
Hubble parameter h and the power 
index n. 

Here, as in any method for es- 
timating the PS, the recovered PS 
is sensitive to the assumed obser- 
vational errors, that enter as well 
the correlation matrix R. To alle- 
viate this problem, we extend the 
method such that also the magnitude 
of these errors is determined by the 
likelihood analysis. This is done by 
adding free parameters that govern 
a global change of the assumed er- 
rors, in addition to modeling the PS, 
and provides some reliability check 
of the magnitude of the errors. We 
find, for both catalogs, a good agree- 
ment with the original error esti- 
mates, thus allowing for a more re- 
liable recovery of the PS. 

3 Results 

Figure |l] shows, as a typical ex- 
ample, the results for the COBE- 
normalized fiat ACDM family of 
models, with a tensor component in 
the initial fluctuations, when setting 
h = 0.6 and varying Q and n. Shown 
are ln£ contours for the SFI cata- 
log and for Mark III. As can be seen 
from the elongated contours, what 
is determined well is not a specific 
point but a high likelihood ridge, 
constraining a degenerate combina- 
tion of the parameters roughly of the 
form fin 3 ' 7 = 0.59 ± 0.08, in this 
case. The corresponding best-fit PS 
for the two catalogs is presented as 
well, with the shaded region illustrat- 
ing the 90% confidence region ob- 
tained from the SFI high-likelihood 
ridge. 

These results are representative 
for all other PS models we tried. 
For each catalog, the different mod- 
els yield similar best-fit PS, falling 
well within each others formal un- 
certainties and agreeing especially 
well on intermediate scales (k ~ 
O.lh Mpc -1 ). The similarity of the 
PS obtained from SFI with that of 
Mark III, which is seen in the fig- 
ure, is illustrative of the other mod- 
els as well. This indicates that the 
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Figure 1. Likelihood analysis results for the COBE-normalized flat KCDM 
model with h = 0.6. Shown are ln£ contours in the f2 — n plane for SFI (top 
left panel) and for Mark III (top right). The best-fit parameters for SFI and 
Mark III are marked, on both, by 'S' and 'M' respectively. The lower panel 
shows the corresponding maximum- likelihood PS for SFI (solid line) and for 
Mark III (dashed). The shaded region is the SFI 90% confidence region. The 
three solid dots mark the PS calculated from Mark III by Kolatt and Dekel 
(1997), together with their quoted la error-bar. 



peculiar velocities of the two cata- 
logs, with their respective error es- 
timates, are consistent with arising 
from the same underlying mass den- 
sity PS. This does not preclude possi- 
ble differences that are not picked up 
by this statistic, but can be viewed as 
another indication of the robustness 
of the results. Note also the agree- 
ment with an independent measure 
of the PS from the Mark III catalog, 
using the smoothed density field re- 
covered by POTENT (the three dots; 
Kolatt & Dekel 1997). 

The robust result, for both cat- 
alogs and all models, is a rela- 
tively high PS, with P(fc)C 1 ' 2 = 
(4.5 ± 2.0) x 10 3 (h^Mpc) 3 at k = 
0.1 h Mpc -1 . An extrapolation to 
smaller scales using the different 



CDM models gives o 8 n 0s = 0.85 ± 
0.2. (Such values are also obtained 
when assuming for the PS the T 
model or the generalized CDM mod- 
els with a free amplitude.) The 
high-likelihood ridge is a feature of 
all COBE-normalized CDM models. 
The general constraints on the com- 
bination of cosmological parameters 
is of the sort H h(, ^ n" = 0.6 ± 0.2, 
where fi = 1.3 and v = 3.7, 2.0 for 
flat ACDM models with and with- 
out tensor fluctuations respectively. 
For open CDM, without tensor fluc- 
tuations, the powers are /i = 0.9 
and v = 1.4. For the span of models 
checked, the PS peak is in the range 
0.02 < k < 0.06/i Mpc -1 . The shape 
parameter of the T model is only 
weakly constrained to T = 0.4±0.2. 
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Figure 2. Constraints in the Q — Ha plane coming from a likelihood analysis 
of the SFI peculiar velocities (the approx. vertical contours) and from the 
high-redshift SNe la (the diagonal contours; Perlmutter etal. 1998). The 68, 
90 and 99% confidence regions are shown for both. The estimated worst-case 
systematic error of the SNe results is denoted by the corresponding dashed 
line. The SFI likelihood contours are for the case of n = 1, h = 0.65. The 
shifted dashed lines illustrate the estimated effect of changing the values of 
these parameters (Zehavi & Dekel 1998). 



These error-bars are crude, reflect- 
ing the 90% formal likelihood uncer- 
tainty for each model, the variance 
among different models and between 
catalogs. Care should also be given to 
possible systematics that could still 
plague the results, arising perhaps 
from non-linear effects or some pecu- 
liarities in the data (Freudling etal. 
1998). 

4 Further Analysis 

We have recently extended the 
analysis of COBE-normalized CDM 
models to models with general values 
of Q and n A (Zehavi & Dekel 1998). 
Although Q A comes in only indi- 
rectly through the COBE normaliza- 



tion, such results are particularly in- 
teresting as they can be compared 
to other constraints in the f2 — f^A 
plane, such as the recent results com- 
ing from high-redshift type la super- 
novae (SNe la; Perlmutter et al. 1998; 
Riess etal. 1998). 

Figure ^ illustrates such con- 
straints in the Q — Qa plane, show- 
ing ln£ contours for the SFI cata- 
log, for fixed values of n = 1 and 
h = 0.65. The peculiar velocity anal- 
ysis appears to constrain an elon- 
gated ridge in this plane of a nearly 
fixed Q and varying Qa- As demon- 
strated in the plot, a change in the 
values of n and h essentially shifts 
the ridge toward either a higher or 
lower Q, for smaller and larger val- 
ues of these parameters, respectively. 
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This is another manifestation of the 
degeneracy between these parame- 
ters mentioned earlier. The accept- 
able range of these parameters is 
therefore needed to be determined by 
other external constraints. The re- 
sults of the same analysis applied to 
the Mark III catalog are here as well 
fairly similar to the SFI ones, except 
for a slightly stronger preference to- 
ward smaller values of Q\ . 

The confidence contours in this 
parameter plane obtained by the Su- 
pernova Cosmology Project (Perl- 
mutter et al. 1398) are sketched as 
well in Figure El (These results are 
consistent with the findings of the 
High-z Supernova Search Team re- 
sults, Riess etal. 1998.) Taking into 
consideration concurrently these two 
independent sets of constraints seem 
to imply a considerable contribution 
from both Q and Qa- While the 
range of models consistent with the 
high-z SNe findings alone includes 
low Q + low fl\ models (and even 
< 0), the peculiar velocities anal- 
ysis appears to rule out these mod- 
els, and taken together the two make 
a stronger case for a positive cosmo- 
logical constant. 

Work in progress includes an at- 
tempt to merge the Mark III and SFI 
catalogs and perform a joint likeli- 
hood analysis, looking specifically at 
the cross-correlations of the two data 
sets, that may entail valuable infor- 
mation. Another aim is to perform 
a similar analysis on other veloc- 
ity data, such as velocities of galaxy 
clusters (e.g, Smith etal., these pro- 
ceedings) or SNe la velocities, which 
probe out to larger scales with rela- 
tively high accuracy. Lastly, an inter- 
esting prospect is to do a simultane- 
ous analysis of velocity data together 
with other kinds of data, like red- 
shift surveys and CMB experiments 
(see Webster etal. 1998; Lahav & 
Bridle, these proceedings). The dis- 
tinct types of data complement one 
another, each constraining different 
combinations of the cosmological pa- 
rameters, and together may remove 
the degeneracy and set tight con- 
straints. 



Acknowledgments 

I thank my collaborators in this 
work L.N. da Costa, A. Dekel, A. 



Eldar, W. Freudling, R. Giovanelli, 
M.P. Haynes, Y. Hoffman, T. Ko- 
latt, J.J. Salzer, G. Wegner and S. 
Zaroubi. This work was supported by 
the DOE and the NASA grant NAG 
5-7092 at Fermilab. 



References 

da Costa L.N., Freudling W., Weg- 
ner G., Giovanelli R., Haynes M.P., 
Salzer J.J., 1996, ApJ, 468, L5 
da Costa L.N., Nusser A., Freudling 
W., Giovanelli R., Haynes M.P., 
Salzer J. J., Wegner G., 1998, MN- 
RAS, 299, 425 

Eldar A. etal, 1998, in preparation 
Freudling W., Zehavi I., da Costa 
L.N., Dekel A., Eldar A., Giovanelli 
R., Haynes M.P., Salzer J.J., Wegner 
G., Zaroubi S., 1998, ApJ, submitted 
Haynes M.P. etal., 1998, in prepara- 
tion 

Kaiser N., 1988, MNRAS, 231, 149 
Kolatt T., Dekel A., 1997, ApJ, 479, 
592 

Kolatt T., Dekel A., Ganon G., 
Willick, J. A., 1996, ApJ, 458, 419 
Perlmutter S. etal., 1998, LBNL 
preprint 41801 

Riess A.G. etal., 1998, AJ, in press 
Strauss, M.S., 1998, in Dekel A., Os- 
triker J. P., eds., Formation of Struc- 
ture in the Universe. Cambridge Uni- 
versity Press, Cambridge, p. 172 
Strauss M.A., Willick J.A., 1995, 
Phys. Rep., 261, 271 
Webster M., Bridle S.L., Hobson 
M.P., Lasenby A.N., Lahav O., 
Rocha G., 1998, ApJL, submitted 
Wegner G. etal., 1998, in prepara- 
tion 

Willick J. A., Courteau S., Faber 
S.M., Burstcin D., Dckcl A., 1995, 
ApJ, 446, 12 

Willick J. A., Courteau S., Faber 
S.M., Burstcin D., Dekel A., Kolatt 
T., 1996, ApJ, 457, 460 
Willick J. A., Courteau S., Faber 
S.M., Burstein D., Dekel A., Strauss 
M.A., 1997, ApJS, 109, 333 
Zaroubi S., Zehavi I., Dekel A., Hoff- 
man Y., Kolatt T., 1997, ApJ, 486, 
21 

Zehavi I., Dekel A., 1998, in prepa- 
ration 



